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, Spin-transfer torques occur in magnetic heterostructures because the transverse component of a 

spin current that flows from a non-magnet into a ferromagnet is absorbed at the interface. We 
demonstrate this fact explicitly using free electron models and first principles electronic structure 
calculations for real material interfaces. Three distinct processes contribute to the absorption: (1) 
spin-dependent reflection and transmission; (2) rotation of reflected and transmitted spins; and 
' (3) spatial precession of spins in the ferromagnet. When summed over all Fermi surface electrons, 

these processes reduce the transverse component of the transmitted and reflected spin currents to 
nearly zero for most systems of interest. Therefore, to a good approximation, the torque on the 
magnetization is proportional to the transverse piece of the incoming spin current. 
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I. INTRODUCTION 



When a current of polarized electrons enters a ferromagnet, there is generally a transfer of angular momentum 
between the propagating electrons and the magnetization of the film. This concept of "spin transfer" was proposed 
independently by SlonczewskiEJ and Bergerel in 1996. Experiments soon followed where anomalies in the current-voltage 
characteristics of magnetic heterostructures were interpreted as evidence for spin transfer .tl Unambiguous confirmation 
came when the phenomenon of giant magnetoresistanceu was used to detect magnetization reversal in ferromagnetic 
multilayers with large current densities flowing perpendicular to the plane of the layers otjl3 Subsequently, spin transfer 
has been implicated to explain the observation of spin precession for high-energy, spin-polarized electrons that traverse 
a magnetic thin filmEj and enhanced Gilbert damping in magnetic multilayers compared to one-component magnetic 
filmsfl More experiments may be expected in the future because spin transfer is expected to play an important role 
in the nascent field of "spin electronics" .E3 

Theoretical work on spin transfer generally falls into one of three categorifffi. |r Qpe group of articles focuses on 



£*"»■ ' deriving and solving classical equations of motion for the magnetizationOO'tiJpliSEj These studies generalize the 
0*\ Landau-Lifshitz equation to take account of spin currents, spin accumulation,!!!] and the mechanical torques which 
necessarily accompany (spin) angular momentum transfer. A se cfrind group of articles generalizes charge transport 
theory to take account of spin currents and spin relaxation.Ea£jc3'EM3 These theories compute the spin-transfer 
torques that serve as input to the magnetization calculations. The torque can be computed phenomenologically, 
or from the Boltzmann, Kubo or Landauer formalisms that incorporate quantum mechanical information explicitly. 
Finally, there are articles that report quantum mechanical calculations of the parameters that serve as input to the 
transport theories. The model studies of SloncaewskiEl and BergerB are of this sort, as are the first-principles, electronic 
structure calculations pf. Xia and co-workers Ea 

In a previous paper the present authors used a 2 x 2 matrix Boltzmann equation to compute spin currents, spin 
accumulation, magnetoresistance, and spin-transfer torques in a Co/Cu/Co multilayer with non-collinear magnetiza- 
tion. The physics of spin transfer entered this semi-classical, kinetic theory calculation through quantum mechanically 
derived matching conditions imposed at each ferromagnet/non-magnet interface. Specifically, we took account of a 
reflection mechanisma that arises because the interface reflection and transmission amplitudes for polarized electrons 
are spin dependent. We also took account of an averaging mechanism! that arises because conduction electron spins 
precess around the magnetization vector in each ferromagnet. The present work was motivated originally by two 
1 assumptions we made to simplify the Co/Cu/Co calculations. First, we set to zero the transverse component of the 
spin of the conduction electron ensemble in each ferromagnet. Second, we disregarded the phase of the reflection and 
transmission amplitudes. As best we can determine, the same assumptions are implicit in the Landauer-type model 
calculations reported in Ref. [l^ and Ref. Therefore, before calculations of this sort are carried very much further, 
it seemed appropriate to look more carefully into the correctness of these assumptions. As we will the spin transfer 
process is more subtle and complex than previously managed. 

In this paper, we analyze quantum mechanically the fate of a polarized current that enters a ferromagnet from 
a metallic non-magnet. Using both the free electron model and first principles electronic structure calculations, we 
conclude that the assumptions in question are largely justified. An important point is that the spin of an electron 
generally rotates when it is reflected or transmitted at an interface. This separates the reflection mechanism into two 
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pieces. A spin-filter effect reduces the transverse spin component of each electron individually. A further reduction 
occurs when we sum over all Fermi surface electrons because substantial phase cancellation occurs when the distri- 
bution of spin rotation angles is broad. As for the mechanism we called "averaging" in Ref.E2l, cancellation occurs 
because electrons have different precession frequencies. This leads to an asymptotic, oscillatory, power-law (rather 
than exponential) decay of the transmitted transverse spin component. Putting everything together, we find that 
(except in very exceptional cases) the transverse spin current is almost completely absorbed within a few lattice 
constants of the interface. None, or very little, is reflected or transmitted. As a result, the spin-transfer torque is very 
nearly proportional to the transverse piece of the incident spin current. 

The plan of this paper is as follows. In Section II, we define the basic variables of spin transport and establish 
our notation. Section III analyzes the spin current and spin-transfer torque near a magnetic/non-magnetic interface 
using a free electron model for both materials. Section IV generalizes the analysis of Section III to the case of real 
materials. We summarize our results in Section V. 



II. BACKGROUND 



To help introduce the theory of spin transport, it is useful first to set down the familiar equations of particle 
transport. These involve the number density, 



and the number current density, 



i(r) = iMr), (1) 



j(r) = Re^C(r)vfe(r), (2) 



where v = — (ih/m)V is the velocity operator and ipi !<T (r) is an occupied single particle wave function with state index 
i and spin index a. The continuity equation, 

dn , . 

V. J + ^0, (3) 

expresses the conservation of particle number. In this paper, we will be interested exclusively in steady-state situations 
where the time derivative in (|J) is zero. Not far from equilibrium, the current takes the phenomenological form, 

j = (<r/e)E - DV5n, (4) 

where Sn = n — n cq is the deviation of the number density from its equilibrium value, E is an electric field, a is the 
conductivity, and D is a diffusion constant. The latter two are second rank tensors in the general case. 
For the spin degree of freedom, the analogs to ([!]) and (Q) are the spin density, 

m ( r ) = ^ M ( r ) ' ( 5 ) 



and the spin current density 



Q(r) = R e (r) s„y ® v </w (r)] , (6) 



where s = (h/2)<r and er is a vector whose Cartesian components are the three Pauli matrices. The spin current is a 
tensor quantity. The left index of Qij(r) is in spin space and the right index is in real space. Spin is not conserved so 
the analog of (||) generally has non-zero terms on the right hand side. For our problem, 

_ — 9m dm 

V • Q + — = + n ext 7 

9t t u 



where n ext is an external torque density, V • Q = dkQik and 6m = (|m| — m cq )rh is the so-called spin accumulation 





The first term on the right side of (0) accounts for the transfer of angular momentum between the spin current and 
the lattice due to spin-flip. This process, with relaxation time T|j_, changes the magnitude of the local spin density, 
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but not its direction. The second term on the right side of (Q) describes all external torques that act to change the 
direction of the local magnetization. For example, the Landau-Lifshitz-Gilbert torque density, 

n ox t = -(gnB/K)m x B cff + am x m (8) 

includes an effective field B e ff and phenomenological damping. The effective field is due to exchange, anisotropies, 
and any external fields that might be present. 

To study magnetization dynamics, we merely rearrange ([?]) to 

dm 

— - = n c + n cxt (9) 
at 

where 

n c = — V Q (10) 

is the current-induced contribution to the torque density. The divergence theorem then shows that, apart from spin- 
flip, the torque on the total magnetization in a volume V arises from the net flux of spin current into and out of the 
surface S that bounds V . Phcnomcnologically, the spin current is driven by drift and diffusion: 

Q lk = afE k - DPdkSrm (11) 

As in (|), we assume the simplest form for the spin transport coefficients. That is, we use the vector <x p = (a^ — a{) m 
rather than a third rank tensor and the scalar D p rather than a fourth rank tensor. The conductivities o~\ and <7j 
refer to majority and minority electrons, respectively. 

In a non-magnet, <7j = <7[ and a spin current arises only if there are regions of the metal where there is a gradient 
in the spin accumulation, Sm(r). This implies that the spin density m(r) and the spin current density Q(r) are only 
indirectly related to each other. They need not be collinear and are generally not proportional. In a ferromagnet, an 
electric field produces a current of polarized spins simply because cx-\ 7^ cr| . This spin current is modified by gradients 
in spin accumulation also. However, the second term on the right side of ( |TT|) is valid (at most) when the direction 
of the ferromagnetic magnetization is uniform in_space. Corrections are necessary when the magnetization rotates 
continuously in space, e.g., inside a domain wall.LH 

With this background, the remainder of this paper is devoted to a detailed analysis of the fate of a spin polarized 
current that flows from a metallic non-magnet into a metallic, single-domain ferromagnet through an ideal, flat 
interface. Specifically, we point the particle current density vector j along positive x, we point the ferromagnetic 
magnetization vector M along positive z, and we fix the interface at x = 0. Figure |] shows three possible steady 
states of pure current polarization in the non-magnet and the associated non-zero component of the spin current 
density tensor. For each case, we let only one component of Q ax be non-zero. Q zx 7^ corresponds to longitudinal 
(parallel to M) current polarization. Q xx ^ or Q yx 7^ correspond to transverse (perpendicular to M) current 
polarization. To produce an "incident" polarized current in the non-magnet, it is sufficient that the current flow into 
the non-magnet from an adjacent ferrmnagnet and that the thickness of the non-magnet be small compared to the 
non-magnet spin-flip diffusion length.EH For this reason, magnetic multilayer structures are the rule in most spin- 
transfer experiments. We refer the reader to Ref. [2^ for some insight into the polarization process for the Co/Cu/Co 
system. 

Figure |l| also indicates that, of the three incident states of pure current polarization shown, only Q zx transmits 
into the bulk of the ferromagnet. The magnet absorbs the transverse components. Furthermore (see below), almost 
none of the transverse spin current reflects from the interface. Therefore, if we choose a rectangular pillbox that 
just straddles the interface, the divergence theorem discussion below (|Io|) implies that a current-induced spin-transfer 
torque is exerted on the interfacial magnetization. To be more precise, Figure ^illustrates such a pillbox and incident, 
reflected, and transmitted charge current density vectors. Integrating the steady state (p = 0) version of the continuity 
equation (||) over the pillbox gives 

= (j in -j tr +r f )-Ax (12) 

where A is the area of the interface. Eq. (|l2|) says that the incoming flux j ln • Ax minus the outgoing flux j tr • Ax + 
j rcf • (—Ax) equals zero. The reflected flux has a minus sign relative to the transmitted flux because it passes through 
the opposing face of the pillbox. 

Ignoring spin flip, the same integration applied to ( |l0| ) yields 

N c = (Q in - Q tr + Q rcf ) ■ Ax ~ Q [ l ■ Ax (13) 
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FIG. 1: Three states of spin current scatter from an interface. The current flows from left-to-right, from the non-magnet into 
the ferromagnet. Q zx is longitudinal (parallel) to the magnetization M. Q xx and Q yx are transverse to M. Only Q zx can be 
non-zero in the bulk of the magnet. The transverse spin currents are absorbed in the interfacial region. 
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FIG. 2: Interfacial pillbox used as the integration volume when the divergence theorem is applied to ^ and ( |lo|) to derive jl^) 
and dl3|). 



where Q in , Q rcf , and Q tr are the spin current density (|^) computed using incident state, reflected state, and trans- 
mitted state wave functions. Eq. says that the incoming spin flux Q m • Ap,minus the outgoing spin flux 
Q tr • Ax + Q rcf • (—Ax) equals the torque on the magnetization inside the pillboxBil The torque N c is a vector in 
spin space because we have contracted the space index of the spin current density with the space vector x. The 
approximate form on the right of Eq. (Jl^) says that the torque is proportional to the tranverse part of Q; n . That is 
the main message of this paper. The following sections are devoted to a demonstration that the tranverse transmitted 
and reflected spin currents do indeed disappear in the immediate vicinity of the interface. 



III. FREE ELECTRONS 



In this section, we compute the spin current near the interface of a non-magnet and a ferromagnet assuming that a 
free electron description is adequate for the conduction electrons in the non-magnet and also for both the majority and 
minority conduction electrons in the ferromagnet. We do this in the interest of analytic simplicity and also because 
some authorstLMLa believe this is a fair representation of reality for the purposes of transport calculations. 

We first work out the problem of one electron scattering from a planar interface to determine the amplitudes for 
reflection and transmission. They turn out to be spin dependent. As first shown by Slonczewskijd this fact alone 
generates a "spin-filter" torque because the wave function for an incident electron with a non-zero component of spin 
transverse to M can always be re-expressed in terms of up and down spin components. 

The actual current polarization in the metal is obtained by summing over the full distribution of conduction 
electrons. This introduces two effects. The first arises because the reflection amplitude for free electron interface 
scattering is complex. This means that the spin of an incoming electron rotates upon reflection. The cancellation 
which occurs when we sum over all these spin vectors reduces the net transverse spin current because reflection and 
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transmission both contribute to the outgoing flux from the interface region. A second effect arises because up and 
down spin electrons on the Fermi surface with the same wave vector in the non-magnet no longer have the same 
wave vector when they transmit into the ferromagnet. The two states are coherent, so precession in space (rather 
than time) occurs. The precession frequency is different for electrons from different portions of the Fermi surface. 
Therefore, when we sum over all conduction electrons, almost complete cancellation of the transverse spin occurs after 
propagation into the ferromagnet by a few lattice constants. 



A. Spin currents for a single electron 



Let us choose the spin quantization axis to be parallel to the magnetization of the ferromagnet. Then, in the 
non-magnet, the wave function for an electron whose spin points in an arbitrary direction can always be written as a 
linear combination of spin up and spin down components. Specifically, 

ip ia = [cos^e-^ 2 |T)+sini0e^/ 2 ||)]e l ' c ^e lq - R (14) 

represents a free electron propagating toward the interface in Figure [j] with its spin pointed in the direction (9, </>) 
with respect to M. We are interested in conduction electrons so the wave vector k = (k x , q) satisfies hk 2 /2m = Ep. 
The spatial variable is r = (x,R). As the notation indicates, ([l4|) is the incident state for a scattering problem th 
determines the wave function for the entire system. The latter describes a steady-state situation like current flow.l 
Like the incident state ( |l4| ) , the complete scattering state can also be written as a linear combination of spin up and 
spin down components: 
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In detail, 

ipl = cosine- 1 */ 2 |T) 
il>l =sini(2e^/ 2 |!) 



V^Vt+V-i- (15) 



{e tk * x + J R T e- lfc ^)e iq - R x < 

T ^ e iklx e iq-R x> q 

(e lk * x + R ie - lk * x )e lcl R x < 

Tj yz4x e iq-R x > o 

(16) 



where R^, R^, Tj, and Tj are the reflection and transmission amplitudes for up and down spin electrons. These 
amplitudes do not depend on the angles 9 and <f>. Notice that the up and down spin components do not propagate 
with the same wave vector for x > 0. The wave vectors differ because their kinetic energy depends on the exchange 



potential energy in the ferromagnet. The common factor of exp(iq • R) in (16) reminds us that scattering from a flat 
interface conserves the wave vector component parallel to the interface. 

The transmission and reflection amplitudes are determined by the magnitude of the potential step at the interface. 
For a constant effective mass, this step height is parameterized by &F,fc F , and /c F < fc P , the Fermi wave vectors 
for, respectively, electrons in the non-magnet, majority electrons in the ferromagnet, and minority electrons in the 
ferromagnet. The usual quantum mechanical matching conditions yield the real transmission amplitudes 



where k x (q) = ^/k\ — q 2 and k x (q) = (fcp) 2 — q 2 . The reflection amplitudes are real or complex depending on the 
magnitude of the parallel wave vector. They are 

Ml) = t^l T MS if 9 2 <(fcp) 2 (18) 



and 



k x {q) + K{q) 



= feral if < 2 ><^ 2 < w > 
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where n x (q) = J q 



Kp 



The associated transmission and reflection probabilities, 

R CT (g) = \R a (q)\ 2 
T . (?) = W| Tcr( g ) |a i 



(20) 



satisfy R" 7 + T a = 1 and are plotted in Figure || for a slice through the free electron Fermi surfaces defined by 
kp/kp = 1.5 and kp/kp = 0.5. For this case, the transmission probability for majority electrons (dashed curve) 
is unity near the zone center and then falls rapidly to zero near fcp. The minority electrons (solid curve) transmit 
similarly except that T^ falls to zero near fcp. 

It is now straightforward to compute and interpret the incident, reflected, and transmitted number current densities 
and spin current densities. We need only (||), (g), and the appropriate piece of the wave function ([l6|). The incident 
current densities are 



3 X 



v., 



v x 
Ti 

h 



-v x smf cos( 



Q'yx = ^Smflsmf 

h 



(21) 



where v x = Tik x /m. The reflected current densities are 

x\ [ cos2 

\v x \ [cos 2 i0|ii T | 2 -sm : 



•rcf 



ref 

xx 



Q 

^y X 



-\v x \ [cos 2 \9\R^\ 2 +sin : 

h 

2 

h 

1 
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:d\RA 2 ] 



^Isin^Re^^e -10 ] 



(22) 



The transmitted current densities are 



°i ZX 



r: x (r) 



fy X (r) 



^cos 2 ^!^^^^ 2 ^!^! 2 

^icos 2 ^^! 2 -^^ 2 ^^! 2 



4 2 



4- 



(23) 



where v% — Hk^/m. Using (|20|), it is easy to check that j™ = j* r — j r x and = Q z r x — Q r z c x ■ The first relation is 
consistent with ( |l2"| ) because there is no accumulation of charge at the interface. Using (|l3|), the second relation tells 
us that there is no torque associated with the transport of longitudinal spin current. However, a similar relationship 
does not hold for the other two components of Q. There is a discontinuity in the transverse spin current when a spin 
scatters from an interface. According to ([k|), this implies that a current-induced torque acts on the magnetization. 
In fact, three distinct mechanisms contribute to the net torque. 

One source of discontinuity and spin-transfer torque is spin filtering. This occurs when the reflection probabilities 
are spin dependent.EJ To see this, note first that the specific superposition of up and down spin components displayed 
in the incident state wave function ( |l4| ) corresponds to a specific transverse component of the spin vector. If i?| = R[ 
and Tf = Tj. , that specific linear combination is preserved in the reflected and transmitted pieces of the scattering state 
and no discontinuity occurs in the spin current. However, if the reflection and transmission amplitudes differ for up 
and down spin components, the up and down spin content of the spatially separated reflected and transmitted states 
differ from one another. This leads unavoidably to different transverse spin components and thus to a discontinuity in 
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the transverse spin current. Given the structure of (£2|) and (f23|), we use the reflection and transmission probabilities 
in the combination vRJRJ + vTTTT as a measure of the ability of spin filtering to provide spin-transfer torque. The 
next-to-top and next-to-bottom panels in Figure |3] display the required information. 

A second source of transverse spin current discontinuity and spin-transfer torque is spin rotation. This occurs when 
the product R*Ri is not positive real. Specifically, ( p2| ) shows that the transverse components of the reflected spin 
current contain a factor 

i? T *i? i = |i?^ i |e lA ^. (24) 

The phase Atp evidently adds directly to the azimuthal angle <j> used to define the spin direction in the incident state 
vector (|l4]). In other words, the reflected spin direction rotates with respect to the incident spin direction. This is 
an entirely quantum mechanical phenomenon for which there is no classical analog. The bottom panel of Figure ^ 
shows that the range of A<j) can be surprisingly large. Indeed, for this choice of Fermi surfaces, the spin direction 
completely reverses when an electron reflects from the interface at near-normal incidence. There is no corresponding 
rotation for transmitted electrons because Tj and T± are positive real (for free electrons). The resulting discontinuity 
in the transverse spin current leads to a spin-transfer torque that is distinct from spin filtering. 

Finally, a glance at ( p3| ) reveals that spin precession is a third source of spin-transfer torque. Note specially 
the spatially-varying phase factors which appear in the transmitted transverse spin currents because k\ ^ k~ in the 
feppmagnet. Their net effect is spatial precession because Q xx and Q yx simply rotate into one another as a function of 
xEB From (0), such a spatial variation of Q contributes a distributed torque density at every point in the ferromagnet. 
The top panel of Figure [}] shows the range of spatial precession "frequencies" 

Ak = k x - tf x (25) 

for the free electron model of that figure. 



B. Spin currents for a distribution of electrons 



The spin currents relevant to experiment reflect the combined effect of all the conduction electrons. In the most 
general description of transport, it is necessary to keep track of the quantum mechanical coherence between all 
electrons in different eigenstates. However, to model the spin-transfer torque experiments reported to datep'B it is 
not necessary to maintain the coherence between states with different Fermi surface wave vectors. It is sufficient to 
use a semi-classical theory that maintains only the coherence between up and down spin states at each k-point on the 
Fermi surface. Accordingly, we define a 2 x 2 electron occupancy distribution matrix 

f(k,r) = [/(k,r) ( /T( J' r) h{ l r) ) C/t(k,r) (26) 

in terms of the scalar occupancy functions for up and down spins and the spinor rotation matrix 

T7(\c v\-( cos(0/2)e^/ 2 - sin(0/2)e-*/ a \ , , 

U[k ' V) ~ { sin(0/2)e^/ 2 cos(0/2)e^/ 2 J' [ ' ' 

We have suppressed the k and r dependence of 8 and 4> for simplicity. ._. 

Elsewhere, we have solved the Boltzmann equation to find f(k,r) for a typical spin-transfer geometryEj For the 
simple scattering problem treated here, the reflected and transmitted distributions are determined entirely by the 
reflection and transmission amplitudes and the incident electron distribution at the interface between the non-magnet 
and the ferromagnet: f(k) = f(x = 0,k). For this distribution, the semi-classical version of the spin current (|^) is 

h f H 3 k 

=2 / (^ Tr [ f ( k )' T ]® v ( k )- ( 28 ) 

•Ux>0 

The restriction v x > limits the integration to the occupied electron states that move toward the interface. Using 
(p6|), ( p7| ) and the cyclic properties of the trace, we get, e.g., 

h r d 3 k 

QL" = j / 7^/p(kK(k) S m0 k cos<£ k , (29) 

u x >0 
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FIG. 3: Slices through a set of free electron Fermi surfaces. The two middle panels show the Fermi surface for the non-magnet 
and the superimposed Fermi surfaces of the majority (dashed) and minority (solid) states of the ferromagnet. The panel just 
above the magnetic Fermi surfaces is the probability for transmission into the ferromagnet for majority (dashed) and minority 
(solid) electrons. The panel just below the non-magnetic Fermi surface is the probability for reflection back into the non-magnet 
for majority (dashed) and minority (solid) electrons. The bottom panel shows the phase in ( |24| ) acquired by an electron because 
its spin rotates upon reflection. The top panel shows the wave vector difference (03) for a transmitted electron. 



where / p (k) = /|(k) — /j,(k) determines the degree of polarization at each point on the Fermi surface. The angles 
£?k and 0k determine the direction of the spin polarization. Electron states in the immediate vicinity of the Fermi 
surface dominate the transport of charge and spin. Therefore, we write 



/ CT (k)^/ (e k )+ 5CT (q) 



<9e k 



(30) 



where f is the equilibrium Fermi-Dirac distribution function and the partial derivative restricts k to the Fermi surface. 
We write g a (q) rather than g a (\s.) because |k| 2 = k 2 + q 2 = fc|. The equilibrium term does not contribute to the spin 
current. Otherwise, we let d 3 k = d 2 qdk x and use J dkxdfo/dek — l/7i|« x (q)| in (|2rj|). The result is 



t^>0 



(2v) 



rffp(q) sin q cos< 



(31) 
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where 



sp(q) = .9T(q) -5i(q)- 

For Q l y X , change cos </> q to sin <^ q in . 

The reflected spin current due to all the conduction electrons is 

d 3 k 



QI ' Cf(r) = 2 7 (^Tr[Rt(k,r)f(k)R(k,r),r]®v rcf (k) 



«x>o 



where 



R(k) 



i? T (k)e 






i2i(k)e 4 



(32) 



(33) 



(34) 



and v rcf (k) is the velocity of a reflected electron with wave vector k. The r-dependent phase factors in ( |34| ) cancel 
out in (|33| ) so, e.g., 



O ret = 

^xx 



47T 



«x>0 



rf 2 g 

(2^ 



■flip(q) sin6» q |i?*(q)E;(q)| Re 



-j(l/>q-A0q) 



(35) 



where A0 q is the relative phase of the reflection amplitude as in (|24j). For Q'£, change Re to Im in fl3§ 



Finally, the total transmitted spin current is 



Q tr (r) 



Ux>0 



d 3 k 
(2^)3 



Tr [T t (k,r)f(k)T(k,r)cr] ® v tr (k) 



(36) 



where 



T(k,r) 



T T (k)e 



ikf. 



and 



v tr (k) 



T x (k)e lk r 
vT(k) +v^(k) 



(37) 



(38) 



In these formulae, the wave vector for incident states, k, transforms to either k^ or k^ when the electron enters the 
ferromagnet. The average transmitted velocity v tr (k) is defined only at values of q where both spins transmit. A 
comparison of, say, 



Qxx( X ) 



1 

47T 



N2 5p(qJ &m y q ~ —■ ^ e 



Vx>0 



|2« s (q) 



Tf(q)T l (q)e' 



(39) 



with ( p3| ) confirms that ( pq ) is correct with the definitions (37) and (|38|). For Q^ x , change Re to Im in (p9|). 

At this point, we must make a specific choice for g p (k) and the polarization of the incident spin current. Let us 
assume the current is completely spin polarized along +x. This fixes #k = 7r/2 and 0k = 0. For the distribution ([32]), 
we begin with the approximate form 



5 P (q) = a + bv x (q). 



(40) 



The two terms account for interface and bulk-effects, respectively. The velocity-dependent bulk term is familiar from 
textbook treatments of electrical conductivitycZI except, from (|Tl|), gradients in spin accumulation (rather than electric 
potential) drive the spin current in the non-magnet. The constant term,is needed because a spin-dependent chemical 
potential difference A/x across an interface also drives a spin current.liH In this paper, we assume that the interface 
resistance is large (large reflection probability) so we use 



,(q) ~ a = Afj,. 



(41) 
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This is the same approximation that is made in Landauer-type transport calculations. 
With these choices, the incident spin current is 



11 r , . i fc| 



^ xx 2(2tt) 2 J ^ 4(2tt) 2 

o 

The reflected spin currents normalized to the incident spin current are 

Qref <~ 



dqqAf,= -— ^A M . (42) 



and 



)rcf 



2 

fep 



dqq|i?:( g )flx(g)|cosA^ q (43) 



fcp 



^^(^(^sinA^. (44) 



i) 



We get Q™* 7^ because, as discussed earlier, many of the spins rotate upon reflection. On the other hand, the 
sinusoidal factors lead to substantial self-cancellation of the integrals (^) and (^) when the range of Atp q is large 
(see bottom panel of Figure ||) £3 In most cases, we find the total transverse reflected spin current to be very small. 
The normalized transmitted spin currents are 
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Based on the behavior of the transverse reflected spin current, we expect ( J45| ) and (|46|) to decay as a function of x 
because the generally wide range of Ak — k\, (q) — kl (q) (see top panel of Figure ||) ought to induce self-cajacellation 
of the integrals. In fact, like a similar integral that appears in the theory of oscillatory exchange couplingjHa we can 
extract the asymptotic form (x — > oo) analytically using a stationary phase approximation. Only small values of q 
contribute in that instance so for, say, the xx component, we find 



J™ %F - -^^^(0)^(0) ,,; K (47) 
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To understand this result, we note (see the top panel of Fig. ||) that the electron states with wave vectors in an interval 
5q near q = (which share the value Ak ~ fcp) play a special role. These states precess together (coherently) with 
spatial frequency kp — kp. Slow dephasing begins only after a distance x where x8q ~ 1. 

The oscillatory, alg«bfaie,decay exhibited by (|4^) contrasts markedly with the assumption of monotonic, exponential 
decay made by others .Hii3E3 Of course, incoherent scattering processes may be expected to superimpose an exponential 
decay on the algebraic decay we find. The solid curves in Fig. || illustrate the behavior of the transmitted spin current 
( fi"5| ) for three free-electron models. The dashed curves show the asymptotic behavior from (^|). The top panel 
corresponds to Fig. ^| where the Fermi sphere of the non-magnet is significantly smaller than the majority sphere 
and significantly larger than the minority sphere. The middle panel is a situation where the Fermi sphere of the 
non-magnet is identical in size to the majority sphere and both are significantly larger than the minority sphere. 
Finally, the lower panel shows results for majority and minority spheres which are, respectively, slightly larger and 
slightly smaller than the Fermi sphere of the non-magnet. This corresponds to the so-called a s-d model" where the 
conduction electrons bands in the ferromagnet are regarded as slightly split by exchange with localized moments. 

The interfacial "spin-filter" makes each solid curve in Fig. |] differ from unity already at x = 0. The filter is most 
effective when the Fermi surface of the non-magnet is poorly matched with one or both of the Fermi surfaces of 
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FIG. 4: Decay of transverse transmitted spin current as a function of distance from the interface for three free electron models. 
In each panel, the solid curve is the exact result ( fl5| ) and the dashed curve is the asymptotic result (^). Top panel: the 
mismatch is very large between the sizes of the magnetic and non-magnetic Fermi surfaces; fcp|/fcF — 1.5 and kpi/kp — 0.5. 
This is the model used in Fig. ^| Middle panel: the Fermi surfaces are identical for the non-magnet and the majority electrons 
in the magnet; /cft/^f = 1-0 and kp^/kF = 0.5. Bottom panel: an s-d-like model where the mismatch is very small between 
the sizes of the magnetic and non-magnetic Fermi surfaces: k-F-\/kp — 1.1 and kpi/k-p = 0.9. 

the ferromagnet. Owing to (|l7]), this is consistent with our earlier discussion where we identified the transmission 
probability condition T*(q) ^ T (q) as a prerequisite to the action of the spin filter. The subsequent decay of each 
curve in Fig. ^ to zero reflects the distribution of spatial precession frequencies as we have indicated. We have repeated 
these calculations assuming that the distribution function g p is proportional to the velocity term in ( [f"o"|) alone rather 
than the constant term in ( ff"o| ) alone. We find no significant changes from the results of Fig. ||. 

We are now ready to use our computed results to find the spin-transfer torque (|l3|) for free electron models. The 
top and bottom panels of Figure || show the transverse spin space vectors N, Q m • x, Q rcf • x, and Q tr ■ x for the Fermi 
surfaces used in the top and bottom panels of Fig. ||. We have suppressed the contraction with x in the spin current 
labels for clarity. In fact, the vectors for Q rcf and Q tr represent these quantities just at the interface. Therefore, 
the reflected piece includes the dephasing effects of differential spin rotation whereas the transmitted piece does not 
include the dephasing effects of differential spin precession. As we have seen, the latter reduces the transmitted spin 
torque to zero not far from the interface. Therefore, we have drawn the torque vector (for a unit area of interface) so 
N = Q ln + Q rcf . The top panel of Fig. || (large Fermi surface mismatch) shows a significant dephasing of the reflected 
spin current. The bottom panel of Fig. a (small Fermi surface mismatch) shows nearly zero reflected spin current. 
The reflected spin current is exactly zero for the model (not shown) used in the middle panel of Fig. ^. These results 
show that, unless the Fermi surface mismatch is very small, the interface effectively absorbs the entire transverse 
component of incident spin current. This abrupt change in angular momentum is the source of current-induced 
spin-transfer torque at the interface between a ferromagnet and a non- magnet. 

The dashed arc labeled Q sf in each panel of Fig. || is a portion of a circle whose center is the "tail" position for all 
three spin current vectors. The radius of this circle, compared to the length of the vector Q ln , gives an indication of 
the magnitude of the spin filter effect. Quantitatively, the circle radius is proportional to 

k k l 
H = A J dqq \R* mi{q) \ + £ fdqq **®*^ |T T (g)n(g)| ■ (48) 
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FIG. 5: Graphical representation of the interfacial torque and transverse spin currents for two free electron models. The 
re-components are horizontal and y-components are vertical. The horizontal arrow is the incident spin current directed along 
the IE-direction. The dashed arc indicates the reduction in spin current due to the "spin-filter" effect. The thick arrow is the 
reflected spin current at ie = 0. The dashed arrow is the transmitted spin current at x = 0. The thin arrow is the final torque 
taking account of the fact that precessional averagingin the ferromagnet drives Q tr — > after a few lattice constants. Panel 
(a) is the large Fermi surface mismatch model of Fig. H. Panel (b) is the s-d model of the bottom panel of Fig. kl 



With this definition, Q sf measures the magnitude of the total outgoing spin current (reflected plus transmitted) 
without taking phase cancellation into account. This scalar measure of the spin filter is truly meaningful only when 
the reflection and transmission amplitudes are both real and positive, which is not the case. Nevertheless, the dashed 
arcs give some insight into the efficacy of the spin filter mechanism for different free-electron Fermi surfaces. 

The foregoing makes clear that free electron models are useful for building intuition about spin currents and spin- 
transfer torque. However, there is no substitute for first-principles calculations if we are interested in specific material 
interfaces. At the very least, such calculations can be used to judge the correctness of approximate constructs such 
as the s-d model. 



IV. REAL INTERFACES 



In this section we repeat the calculations of Section [II B for ten lattice-matched interfaces between a non-magnet 



and a ferromagnet using a more realistic model of the electronic structure for both. Specifically, we calculate the 
transmission and reflection amplitudes using a linearized- augmented-plane- wave implementation of the local-spin- 
density approximation. The details can be found in Rcf. |3y and Ref. [H| Compared to that earlier work, the 
calculations reported here use a mesh in reciprocal space that is a factor of two denser in each direction. For one case 
(Co/Cu), we checked that no changes in relative spin currents greater than 10 -3 occurred when the mesh was made 
another 2x2 denser. Evanescent states (which decay exponentially away from the interface) play a crucial role in the 
calculation of the reflection and transmission amplitudes. We have ignored them in our spin currents computations. 
Their effect is to change the wave functions in the immediate vicinity (a few atomic layers) of the interface in such 
a way that there is no true discontinuity in the transverse spin current at the interface. As a practical matter, this 
means only that the "interfacial" torque we compute is-in reality-spread out over a few atomic spacings. 

The two middle panels of Fig. ^ show a slice through the Fermi surface of copper and the same slice through the 
majority (dashed lines) and minority (solid) Fermi surfaces of cobalt for the Co/Cu(lll) system. The Fermi surface 
topologies here are much more complicated than the corresponding free electron topologies (cf. Fig. ^). Moreover, 
as the Co minority Fermi surface shows, there can be more than one pair of states for each parallel wave vector. 
Consequently, we supplement every integral over parallel wave vectors with a sum over all possible states that move 
toward the interface for each parallel wave vector. We index these states by n, refer to them as associated with the 
n th sheet of the Fermi surface, and adopt the notation k„ CT = {q, k^} to label Fermi surface wave functions. We drop 
the spin index a in the non-magnet. 

The transverse pieces of the incident spin current for a real interface are 

= 47 J (^E5p( k n) sin0 ( k «)cos^(k„) (49) 
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FIG. 6: Same as Fig. ^| for a real material interface: Co/Cu(lll) with parallel wave vectors in the [110] direction. The two 
middle panels show the Fermi surface for the non-magnet and the superimposed Fermi surfaces of the majority (dashed) and 
minority (solid) states of the ferromagnet. The panel just above the magnetic Fermi surfaces is the probability for transmission 
into the ferromagnet for majority (dashed) and minority (solid) electrons. The panel just below the non-magnetic Fermi surface 
is the probability for reflection back into the non-magnet for majority (dashed) and minority (solid) electrons. The bottom 
panel shows the phase in (JsHJ) acquired by an electron because its spin rotates upon reflection. The top panel shows the wave 
vector difference (|62h for a transmitted electron. 



and 

Q'y X = ^ J ^^5p(k„)sin0(k„)sin0(k n ). (50) 

These differ from the corresponding free electron formulae by the sum over the sheet index n. That sum is restricted 
to the sheets of the Fermi surface where the electrons move toward the interface. 

As before, the efficacy of the spin filter can be judged from the interface transmission and reflection probabilities. 
These state-to-state [n — > n') quantities are 

V n' 2 
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and 
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(52) 



The absolute value is needed in (p2| ) because v n > < and R^ n / must be non-negative. Fig. |] shows the transmission 
and reflection probabilities for one slice through the Co/Cu(lll) Fermi surfaces. 
The transverse components of the reflected spin current are 
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(54) 



Here, the sum over n' is restricted to the sheets of the Fermi surface where the electrons move away from the interface. 
Similar to the free electron case, the dephasing of the reflected transverse spin current is determined by reflection 
phases A0^ n , (q) where 
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(55) 



The bottom panel of Fig. ^ shows that the Co/Cu(lll) phases are both more complicated and exhibit greater 
dispersion than the corresponding free electron results plotted in Fig. ||. 

The transverse_Qieces of Q m and Q rcf written above are. closely related to the mixing conductance, G m i x , introduced 
by Brataas et alS3 and computed recently by Xia et alEB In our notation, 
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This formula is relevant to situations where g p (k„), #(k„), and <fi(k n ) in (|49])-(54) are all constants-a restriction 
implicit in the Landauer description of transport. In that case, the real and imaginary parts of G mlx are proportional 
to the xx and yx components of Q m + Q rcf . From (|T3|), the, latter is proportional to the spin-transfer torque if we 
neglect the transverse part of the transmitted spin current P^H 22 ! For the systems treated by both of us, our numerical 
results for spin-transfer torque agree semi-quantitatively with the mixing conductance calculations of Xia at al. 
The transverse transmitted spin currents are 
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Apart from the sums over n' and n" (both restricted to sheets of the ferromagnetic Fermi surfaces where electrons 
move away from the interface), these formula are less simple than the corresponding free electron results ( ]45| ) and 
( [i"6| ) for two reasons. First, the transmission amplitudes, T ann i, are complex rather than real. Second, the Bloch wave 
functions ?/v(R, x, k ncr ) have a non-trivial dependence on the spatial variable R parallel to the interface plane. For 
the latter reason, the transmitted spin currents contain a factor $ n //„/(q, x) defined by 
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This yields 
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when free electron wave functions are used in (J59|). Otherwise, $ n " n '(q, x) is a complex, periodic function of x with 
period equal to one layer spacing. Thus, it can be calculated once and propagated from layer to layer. A related 
factor enters the reflected spin currents (|57|) and J58|). However, because the spin up and spin down wave functions 
are the same in the non- magnet, it reduces to the velocity factor in the numerator of those expressions. 
Given the foregoing, it is sensible to define transmission phases A<^ n , (q) so that 

Tf nn ,T ;nn ,$ n , w (q, x = 0) = | Tf, m ,T w (q, x = 0)| e lA ^>^\ (61) 

This tells us that, unlike free electrons, the spins of Bloch electrons generally rotate when they transmit through a 
real material interface. If the distribution of transmission phases is broad, substantial cancellation of the transmitted 
spin current occurs at x — when we sum over all transmitted electrons. This effect is independent of the spin filter, 
which also acts at x = 0. 

Any spin current that survives to propagate into the ferromagnet rapidly disappears due to differential spatial 
precession. The (generalized) spatial precession frequency is determined by the difference in wave vector for different 
sheets of the majority and minority Fermi surfaces: 

Ak n i n ii = fe^/j — ■ (62) 

The top panel of Fig. ^ illustrates the distribution of Afc„'„" for a Fermi surface slice of Co/Cu(lll). The large 
dispersion seen there suggests that the spin current decays very quickly in the ferromagnet. This is confirmed by 
Fig. [?], which shows the computed decay of the transverse spin current for three interface orientations of Co/Cu. The 
non-zero value of the dashed curves at x = shows that a large amount of rotation occurs upon transmission. The 
Fermi surfaces are more complicated than the free electron models, so the initial decay is more complicated also. 
Nevertheless, both the (111) and (110) orientations settle into behavior that is readily characterized as a damped 
precession. The amplitude of the precession for the (100) orientation is so small that it is difficult to see whether it is 
precessing or not. In general, there could be several decaying precessions with different precession rates and different 
amplitudes. 

It is worth noting that none of these curves (or analogous curves for the other material pairs we have studied) 
resembles the the bottom panel of Fig. || appropriate to the s-d model. This lack of agreement is consistent with the 
fact that essentially all the Fermi surface wave functions in third-row ferromagnets contain as much "localized" 3d 
character as "delocalized" 4s character. c3 

Figure^ graphically summarizes our first-principles spin current calculations for ten different interface combinations. 
The vectors labeled Q ref and Q tr correspond to x — and reflect the effect of spin filtering and spin rotation only. 
Q rof is very small and, as we have emphasized, Q tr — * after a few lattice constants. Therefore, the torque per unit 
area of interface is Q m + Q rcf ~ Q m . Due to spin filtering, differential spin rotation, and differential precession, nearly 
all of the incident transverse spin current is absorbed in the immediate vicinity of the interface. For Co/Cu, Fe/Ag, 
and Fe/Au, the spin-filter accounts for somewhat more than half of the effect and the interface dephasing for the rest. 
For Ni/Cu and Fe/Cr the spin filter effect is weaker. For Ni/Cu, the decay of the precessing transmitted spin current 
plays a large role. 

Of course, our calculations pertain to ideal, lattice-matched interfaces. A variety of effects make the interfacial 
absorption of transverse spin even more efficient. We have mentioned already that scattering in the ferromagnetic 
layer increases the rate of decay of the precession. Steps at the interface lead to increased dephasing for both reflection 
and transmission. For thin layers where the decay of the precession might not be complete, the dephasing on passing 
through the second interface generally leads to a further decay of the transverse spin current. Thickness fluctuations 
further reduce the spin current. 

The interface, torque we compute is interesting because of the recent demonstrations of current-induced magnetiza- 
tion switching .ErtJ However, the material pair that optimizes the switching is not determined by the conversion of the 
spin current into a torque. This process is the same for all of the interfaces considered. For the ideal interfaces consid- 
ered here, the optimum choice depends on the ability of the material pair to generate a spin current in the first place. 
This depends on the spin-dependent interface resistances and the spin-dependept bulk conductivities. The Fe/Au 
and Fe/Ag pairs have the strongest spin dependence of the interface resistance£3 However, in reality the optimum 
combination will likely depend on growth considerations. The general mismatch between the body-centered-cubic Fe 
lattice and the face-centered-cubic Au or Ag lattice will probably lead to poor growth, unless the interface is forced 
to be (100) (where the rotated lattices match quite well) and the number of steps at the interface is kept quite small. 

V. SUMMARY 

In this paper, we used free electron models and first principles electronic structure calculations to study the spin- 
transfer torque that occurs when a spin-polarized current flows from a non-magnet into a ferromagnet through a 
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FIG. 7: Decay of transverse transmitted spin current as a function of distance from the interface for three orientations of 
Co/Cu. For a unit incident transverse polarization, the solid curve in each panel is Q X x(x). The dashed curve in each panel is 

perfect interface. The origin of the torque is a transfer of spin angular momentum from the conduction electrons to 
the magnetization of the ferromagnet. The origin of the angular momentum transfer is the absorption of transverse 
spin current by the interface. We identified three distinct processes that contribute to the absorption: (1) spin- 
dependent reflection and transmission; (2) rotation of reflected and transmitted spins; and (3) spatial precession 
of spins in the ferromagnet. When summed over all Fermi surface electrons, these processes reduce the transverse 
component of the transmitted and reflected spin currents to nearly zero for most systems of interest. Therefore, to 
a good approximation, the torque on the magnetization is proportional to the transverse piece of the incoming spin 
current. 

To be more quantitative, we used the analogy between charge current and spin current to show that a spin current 
flowing in the +x direction (perpendicular to the interface) delivers a torque per unit area 

^ = (Q in - Q tr + Q rcf ) • x (63) 

to a microscopically small region around the interface. Here, Q m , Q tr , and Q rcf are the incident, transmitted, and 
reflected spin currents computed using incident, transmitted, and reflected wave functions. We found the latter 
by solving the one-electron stationary-state scattering problem. In the quasi-classical approximation, the total spin 
current is the sum of contributions from every conduction electron. 

Quite generally, the component of N c parallel to the ferromagnetic magnetization is zero. This is consistent with 
our classical intuition. On the other hand, we found that the transverse components of Q tr and Q rcf relevant to (|63] ) 
are also zero (or nearly so), except in very exceptional cases. This means that the entire transverse spin current is 
absorbed (transferred to the magnetization) in the immediate vicinity of the interface. As indicated above, this is so 
due to spin filtering, differential spin reflection, and differential spin precession. 

The spin-filter effect occurs because the wave function for an incident electron with a non-zero spin component 
transverse to the magnetization spin can always be written as a linear combination of spin-up and spin-down compo- 
nents. Then, because the reflection and transmission amplitudes differ for up and down spins, the up and down spin 
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FIG. 8: Graphical representation of the interfacial torque and transverse spin currents for a series of real interfaces. The 
re-components are horizontal and y-components are vertical. The horizontal arrow is the incident spin current directed along 
the IE-direction. The dashed arc indicates the reduction in spin current due to the "spin-filter" effect. The thick arrow is the 
reflected spin current at ie = 0. The dashed arrow is the transmitted spin current at x = 0. The thin arrow is the final torque 
taking account of the fact that precessional averaging in the ferromagnet drives Q 1 — » after a few lattice constants. 



Cu/Co(lll) 



content of the reflected and transmitted wave functions (which are spatially separated) differ from both each another 
and from the incident state. The spin currents directly encode this information. As a result, the right side of (^) is 
non-zero. This is a one-electron effect that operates independently for every electron. 

The two other effects that help drive the transverse parts of Q tr and Q rcf to zero occur when we sum over the 
entire ensemble of conduction electrons. The first arises because the spin of an electron generally rotates when it 
is reflected or transmitted at the interface between a non- magnet and a ferromagnet. The rotation is non-classical 
and the amount of rotation differs considerably for electrons with wave vectors from different portions of the Fermi 
surface. Phase cancellation occurs when we sum over all electrons. In the end, we find that very little remains of the 
reflected transverse spin current. The cancellation of the transmitted spin current is less dramatic. 

Finally, due to exchange splitting, the electrons that transmit into the ferromagnet possess spin-up and spin-down 
components with the same total energy, Ep, but different kinetic energy and so different wave vectors. This implies 
that each electron spin precesses (in space) as it propagates away from the interface. However, like the spin rotation 
angles, the spatial precession frequency varies considerably over the Fermi surface. Consequently, rapid dephasing 
of the transverse spin components of the individual electrons occurs as the conduction electron ensemble propagates 
into the ferromagnet. The net result is a precessing spin current that damps out algebraically within a few lattice 
constants of the interface. 

Our first principles calculations show that the relative importance of these three mechanisms differs for different 
materials pairs and also for different crystallographic orientations for the same material pair. Nevertheless, the final 
result is the same in all cases: the transverse spin current essentially disappears at the interface. The concomitant 
transfer of angular momentum delivers a torque to the magnetization in the immediate vicinity of the interface. 
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Eq. (13) is equivalent to Eq. (10) in [22]. That relation was written in terms of the magnitude of the current. Thus, the term 
for the reflected spin flux appears there with no minus relative to the transmitted spin flux because the reflected current is 
in the opposite direction to the transmitted current. 

We exploit the fact that the cross section for a time-dependent wave-packet scattering process is determined by the scattering 
amplitude (here, the transmission and reflection amplitudes) of the corresponding stationary-state, wave-function matching 
problem. See, e.g., K. Gottfried, Quantum Mechanics (Addison- Wesley, Reading, MA, 1966), Section 12. 
Spatial precession occurs when a coherent superposition of up and down spin states have the same energy (Ef in this case) 
but different wave vectors. The more familiar case of temporal precession occurs when a coherent superposition of up and 
down spin states have the same wave vector but different energies. In both cases, an accumulation of phase changes the 
orientation of the spin vector. 
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